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ABSTRACT 


This  report  presents  the  stability  and  control  derivatives 
of  a low  speed  Navy  trainer  aircraft,  the  T-2B,  as  determined 
from  flight  test  data.  The  maximum  likelihood  method  is  used 
for  the  estimation  of  the  aerodynamic  derivatives,  instrument 
biases  and  random  errors,  and  gust  effects.  The  two  flight 
conditions  considered  are:  (a)  low  speed,  gear  down,  flaps 

deflected,  sea  level,  and  (b)  high  speed,  gear  and  flans  up, 
10,000  feet  altitude.  It  is  shown  that  it  is  necessary  to 
include  control  input  bias  and  in  some  cases  gust  effects  to 
get  reliable  results.  Techniques  for  handling  input  bias  and 
random  noise  in  input  measurement  are  given. 

The  identified  models  in  each  case  are  verified  through 
several  prediction  runs. 
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I.  INTRODUCTION  AND  BACKGROUND 


The  Naval  Air  Development  Center  (NADC)  has  conducted  ex- 
tensive flight  tests  of  a T-2B  aircraft.  The  object  of  these 
tests  is  to  provide  a comprehensive  data  base  of  measurements  of 
the  aircraft  response  which  could  be  used  to  estimate  stability 
and  control  derivatives.  The  data  has  been  supplied  to 
Systems  Control,  Inc.  (SCI)  for  processing  with  an  advanced 
maximum  likelihood  parameter  identification  program.  This  report 
details  the  results  of  that  processing. 

The  T-2B  aircraft  used  for  the  flight  tests  was  extensively 
instrumented,  including  measurements  of  all  attitude  angles, 
rates,  and  accelerations.  The  aircraft  responses  were  excited 
by  several  types  of  inputs  at  different  flight  conditions  in 
order  to  isolate  the  input  effectiveness  in  allowing  accurate 
derivative  estimates.  Specifically,  the  flight  data  records 
supplied  by  NADC  are  at  two  flight  conditions  and  have  excitation 
of  both  short  period  and  phugoid  modes  in  the  longitudinal  motion. 
Flight  condition  1 is  an  approach  configuration  at  236  ft  sec  1 
and  sea  level  with  gear  down  and  flaps  deflected.  Flight  condition  2 
is  at  679  ft  sec  1 and  10,000  ft  altitude  in  clean  configuration 
[see  Table  1],  The  wind  tunnel  values  of  stability  and  control 
derivatives  at  these  flight  conditions  arc  given  in  Table  2. 

To  process  this  data,  a maximum  likelihood  algorithm  was 
used  which  could  identify  random  input  noise  effects  on  the 
aircraft  states.  The  algorithm  used  to  achieve  this  is  highly 
efficient.  Parameters  which  are  not  identifiable  are  isolated 
automatically  during  the  computation.  The  measurement  bias 
and  random  white  noise  errors  are  determined.  The  process 
noise  power  spectral  density  is  identified  in  certain  cases.  The 
results  are  compared  to  wind  tunnel  values  and  to  results  ob- 
tained by  using  simple  output  error  programs,  which  do  not  con- 
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Table  1 

T-2  Flight  Conditions  and  Configurations 


Flight  Condition  1 

1 (Approach 
Configuration) 

mmsmm 

fWiTiRmWiffn 

Mach  No. 

0.212 

0.63 

True  Airspeed  (ft/ sec) 

236 

679 

Altitude  (ft) 

sea  level 

10,000 

Gear  Position 

down 

up 

Flap  Position  (deg) 

16 

0 

Speed  Brake  Position 

closed 

closed 

Weight  (lbs) 

11,000 

11,000 

CG  Position  (t  c) 

20 

20 

Iy  (Slug  - ft2) 

14,600 

14,600 

Ix  (slug  • ft2) 

9,000 

9,000 

Iz  (slug  - ft2) 

19,000 

19,000 

a Trim  (deg) 

4.7 

1.2 

Table  2 

T-2  Dimensional  Longitudinal  Derivative  Estimates 

From  Wind  Tunnel 


Flight  Condi t ion 


1 (Approach 
Gnnfigurat ion) 

2 (Clean 
Configurnt  ion) 

Xu 

(1/scc) 

-0.046 

-0.016 

Xa 

(ft/(sec2  - rad)) 

-23.2 

-36.91 

Zu 

(1/ft) 

-0.00114 

-0.00018 

Zo 

(1/scc  - rad)) 

-0.974 

-2.458 

z6t. 

(l/(scc  - rad)) 

-0.102 

-0.178 

M 

u 

(l/(ft  - sec)) 

0.006 

0.0017 

Ma 

(l/(scc2  - rad)) 

-4.  59 

-22.50 

Ma 

(l/(scc  - rad)) 

-0.53 

-1.57 

Mq 

(i/(scc  - rad)) 

-1.42 

-3.39 

M*c 

(l/(sec2  - rad)) 

K7 

O 

i 

-52.12 

sider  identif iabi lity  problems,  and  with  results  ignoring  process 
noise  (e.g.,  wind  gust  effects  and  noise  in  the  measurement  of 
input) . 

The  comparison  of  various  inputs  shows  clearly  that  certain 
inputs  are  more  efficient  than  others  in  that  they  provide  better 
estimates  of  parameters.  In  general,  the  inputs  of  larger  ampli- 
tude and  lasting  for  longer  periods  of  time  are  desirable.  Also, 
for  a fixed  length  of  the  flight  test  and  input  amplitude,  the 
specially  designed  inputs  are  better  than  conventional  inputs 
such  as  pulses  and  doublets.  It  is  shown  that  it  is  important 
to  estimate  the  bias  in  the  inputs  to  get  accurate  parameter 
estimates . 

The  parameter  estimates  from  different  runs  at  the  same 
flight  condition  show  comparable  mean  values  of  the  estimates  over 
various  inputs.  A prediction  test  (predicting  one  response  from 
the  parameters  estimates  of  a previous  ran)  is  also  shown  to  be 
success  ful . 


Ii.  EQUATIONS  OF  MOTION  AND  MEASUREMENTS 


* < 


In  the  linear  uncoupled  aircraft  model,  the  longitudinal 
equations  of  motion  of  an  aircraft  can  be  described  by  four  states. 
They  are  angle-of -attack  a,  change  in  forward  speed  u,  pitch  rate 
q and  pitch  angle  0.  The  input  is  the  elevator  deflection  6fi. 
The  equations  of  motion  governing  these  states  are  [1] 
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where  ag  is  the  change  in  effective  angle-of -attack  due  to  wind 
gusts.  The  wind  disturbance  is  modeled  as  a first  order  exponen- 
tially correlated  process  (Dryden  Model). 


a * - uia  + w (2) 

g 8 

ui  is  break  frequency  and  w is  a white  noise  process.  We  can  com- 
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There  are  12  parameters  in  the  longitudinal  equations  of  motion. 

They  are  Za,  Zu,  , Xa,  Xu,  Xq,  Ma , Mu,  M^,  Z6e*  X$e,  and  M6g. 

There  are  two  additional  unknown  parameters  associated  with  the  wind 
model-power  spectral  density  of  w and  break  frequency  id.  This  fifth 
order  system  has  five  poles.  There  is  a slow,  lightly  damped  com- 
plex pole  pair  for  the  phugoid  motion  and  another  fast,  highly 
damped,  complex  pair  for  the  short  period  motion.  The  fifth  pole 
(e.g.,  pole  associated  with  the  wind  model)  is  at  w.  If  there  is 
only  short  period  excitation,  the  velocity  is  almost  constant  and 
we  can  use  the  following  four  state  model 
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When  there  is  no  measurement  of  pitch  angle  0 and  the  gust  can 
be  assumed  a white  noise  process,  the  usual  two  state  approximation 
to  the  short  period  motions  can  be  applied. 
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Measurements 


There  are  ten  measurements  of  inputs,  angles,  rates  and 
accelerations  (excluding  speed  brake,  which  is  kept  fixed)  at  a 
sampling  interval  of  0.05  sec.  Table  3 shows  the  measurements  and 
the  functional  relations  to  inputs  and  states.  Five  longitudinal 
measurements  were  chosen  for  identification  of  the  aircraft  deriva- 
tives, unless  otherwise  mentioned.  They  are  angle-of -attack , pitch 
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angle,  pitch  rate,  airspeed  and  normal  acceleration  in  addition  to 
the  elevator  deflection  input.  The  angle-of-attack  is  measured 
using  an  u-vane.  Since  the  a-vane  is  not  located  at  the  C.G.  of 
the  aircraft,  it  also  detects  a component  of  pitch  rate.  The  equa- 
tions governing  the  five  measurments  are 
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V b6’  baz  and  fivc 

white  noise  covariances.  These  parameters  are  included  in  the  set 
to  be  identified. 
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III.  THE  IDENTIFICATION  PROCEDURE 


3.1  INTRODUCTION 

Aircraft  identification  is  the  process  of  determining  aircraft 
stability  and  control  derivatives  from  flight  test  data.  The 
integrated  identification  procedure  is  shown  in  Figure  1.  The 
control  inputs  are  chosen  to  excite  the  desired  modes  of  the  air- 
craft. The  flight  test  is  performed  and  the  measured  output  vari- 
ables are  recorded.  Then  the  identification  algorithm  is  used  to 
process  the  data  and  provide  estimates  of  the  desired  stability 
and  control  derivatives  and  other  unknown  parameters  in  the  instru- 
ment models  or  gust  models. 

Figure  2 indicates  the  various  steps  in  the  data  processing 
stage.  The  input-output  data,  together  with  the  physical  know- 
ledge of  the  system,  are  used  to  formulate  a model  and  the  set  of 
parameters  to  be  identified.  Then,  the  response  data  is  used 
again  to  estimate  the  parameters  in  the  selected  model.  Here, 
one  of  various  techniques  can  be  used.  The  parameter  estimation 
error  covariances  are  also  determined  in  this  stage.  Next,  the 
formulated  model  with  the  identified  parameters  is  verified.  This 
is  the  most  important  stage  in  the  identification  and  is  often 
overlooked.  If  the  model  verification  test  fails,  the  model  may 
have  to  be  reformulated  or  the  flight  test  repeated. 

3.2  CHOICE  OF  PARAMETER  IDENTIFICATION  TECHNIQUE 

Two  identification  techniques  have  been  used  to  identify 
the  stability  and  control  derivatives  as  well  as  bias  and  noise 
covariances  in  instruments  from  flight  test  data. 
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Figure  2 Three-Step  Procedure  lor  Flight  Data  Analysis 


3.2.1  Instrumental  Variables  Approach 


< 


This  is  an  approach  in  which  the  parameters  are  determined 
from  a single  pass  through  the  data.  This  is  a least  squares 
type  approach  except  that  the  bias  of  the  least  squares  estimates 
(because  of  noise  in  the  measurements  of  response  variables) 
is  eliminated  by  suitably  chosen  data  transformations.  The 
technique  is  very  fast,  does  not  require  good  parameter  estimates 
to  ensure  convergence,  and  gives  reasonable  parameter  estimates. 
In  this  application,  the  instrumental  variables  approach  is  used 
together  with  the  wind  tunnel  values  for  start  up  for  the  more 
advanced  and  more  accurate  maximum  likelihood  method.  The  tech- 
nique is  described  in  Reference  2. 

3.2.2  Maximum  Likelihood  Identification 

In  the  maximum  likelihood  approach,  the  likelihood  of  the 
parameters  given  the  measurements  is  maximized.  Conceptual ly , 
this  technique  can  be  summarized  as  follows: 

Find  the  probability  density  functions  of  the  observa- 
tions for  all  possible  comoinations  of  unknown  parameter 
values.  Select  the  density  function  whose  value  is 
highest  among  all  density  functions  at  the  measured 
values  of  the  observations.  The  corresponding  parameter 
values  are  the  maximum  likelihood  estimates. 

The  likelihood  is  a complicated  nonlinear  function  of  the 
mesaurements  and  the  parameters.  Several  techniques  have  been 
used  to  maximize  this  function.  There  are  some  advantages  to 
using  the  Gauss-Newton  gradient  approach  in  this  case  [3].  A 
Kalman  filter  is  used  for  state  estimation  when  process  noise  is 
present . 


i i 


r • 


One  procedure  to  implement  the  Gauss-Newton  optimization 
technique  to  maximize  the  likelihood  function  for  a system 
governed  by  linear  equations  of  motion  is  shown  in  Figure  3. 

The  important  thing  to  note  is  that  the  first  and  second  gradi- 
ents of  the  likelihood  function  can  be  determined  in  terms  of 
the  first  gradient  of  the  innovation  process.  As  has  been  shown 
by  Gupta-Mehra  (3|,  the  innovations  and  their  gradients  with 
respect  to  parameters  can  he  computed  very  efficiently  through 
sensitivity  function  reductions  technique. 

If  there  is  no  process  noise,  the  likelihood  function  re- 
duces to  the  output  error  criterion.  The  details  of  this  method 
are  given  in  Appendix  A. 

3.3  MODIiL  VF.RI  F I CATION 

In  the  model  verification  stage,  we  perform  one  or  more 
tests  to  verify  the  formulated  model  and  the  identified  para- 
meters . 

Several  checks  arc  used  when  the  parameters  are  being  identi- 
fied. The  time  history  of  the  aircraft  responses  computed  from 
the  identified  model  should  be  close  to  the  measured  responses. 

The  difference  (i.e.,  the  innovations)  should  be  a white  noise 
process  and  the  root -mean- square  value  of  this  process  should  be 
close  to  the  noise  in  the  instruments  (slightly  higher  if  there 
is  process  noise).  Also,  the  estimates  of  stability  and  control 
derivatives  obtained  from  different  runs  at  the  same  flight 
condition  should  be  close  to  each  other.  It  is  important  to  note 
the  estimated  covariance  of  the  parameter  identification  errors. 

If  the  standard  deviations  of  estimation  errors  is  of  the  same 
order  as  parameter  values,  then  the  identified  parameters  are 
not  use ful . 


1 1 


FLIGHT  TEST  DATA.  WIND  TUNNEL  VALUES  OF 


AERODYNAMIC  PARAMETERS 


I 

I 


« 


Figure  3 I:low  Chart  of  Maximum  Likelihood  Identification  Program 


All  these  arc,  of  course,  only  the  necessary  conditions  for 
a realistic  model.  It  has  been  suggested  by  several  researchers 
(e.g.,  Akaike  [4|  and  Taylor  [5]),  that  one  of  the  useful  tests 
is  the  prediction  capability  of  the  identified  model.  A schematic 
chart  of  how  this  is  done  to  verify  the  identified  model  of  the 
T-2  aircraft  is  shown  in  Figure  4.  The  measured  control  input 
and  aircraft  response  is  used  to  identify  the  parameter  of  the 
model.  Then  the  aircraft  response  to  another  flight  control 
input,  at  the  identical  flight  condition,  is  measured.  At  the 
same  time,  the  identified  model  is  driven  by  the  same  control 
input.  The  measured  response  of  the  aircraft  will  be  close  to  the 
predicted  response  if  the  identified  model  is  correct. 

3.4  CONSIDERATIONS  IN  THE  USE  OF  THE  MAXIMUM  LIKELIHOOD 

TECHNIQUE 

The  maximum  likelihood  formulation  of  Appendix  A holds  strict 
ly  for  the  case  when  the  "true"  model  is  linear  and  the  measure- 
ment and  process  noise  sources  are  white*  In  actual  engineering 
applications,  this  is  never  the  case.  Therefore,  the  use  of 
the  maximum  likelihood  technique  in  parameter  estimation  problems 
requires  considerable  engineering  judgment.  Many  failures  of  the 
maximum  likelihood  method  arc  directly  attributed  to  a lack  of 
proper  implementation.  Several  considerations  are  important  in 
the  identification  of  aircraft  stability  and  control  coefficients 
using  the  maximum  likelihood  approach: 

1.  Wind  Gust  Effects 


In  many  implementations,  the  wind  gusts  are  not  included 
in  the  model.  This  produces  biased  parameter  estimates. 

In  simpler  terms,  there  is  a consistent  error  in  the  esti- 
mated value  of  parameters.  This  effect  is  included  as 


*The  model  we  refer  to  here  is  the  aircraft  plus  wind  model.  The  g 
acting  on  the  aircraft  usually  cannot  be  approximated  as  white  noi 
A good  approx i mat  ion  of  the  gu  Is  is  white  noise  passed  through  a 
1 i nea r filter. 
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Figure  4 Model  Verification  Via  Independent  Prediction 


detailed  ir  Chapter  II  of  this  report.  By  including  and 
identifying  the  gust  statistical  model,  the  bias  in  the 
parameters  is  reduced. 

2 . Bias  and  Noise  in  Control  Input 

I 

The  elevator  deflection  input  is  usually  measured  at 
the  aerodynamic  surface  using  a potentiometer.  Therefore,  the 
measurement  of  the  input  is  contaminated  by  calibration  biases 
and  random  noise  in  the  potentiometer.  There  arc  additional 
bias  effects  because  of  uncertainty  in  the  trim  condition  of 
the  elevator  deflection.  Notice  that  this  is  an  important 
effect  because  we  take  deviations  from  the  trim  condition  to 
identify  parameters  in  the  aircraft  model.  The  modification  to 
the  maximum  likelihood  method  when  the  inputs  have  bias  and 
random  noise  is  given  in  Appendix  B. 

| 

3 . Identifiabil ity  Problems 

The  maximum  likelihood  method  gives  the  parameter  esti- 
mates as  well  as  the  predicted  error  covariance  of  the  esti- 
mation error.  As  discussed  earlier,  these  error  covariances 
play  an  important  role  in  the  integrated  parameter  identifi- 
cation process.  Special  attention  should  be  given  to  these 
covariances  because  a large  covariance  implies  an 
inaccurate  estimate.  Parameter  estimates  can  be  improved 
by  using  a proper  control  input  in  a subsequent  flight  test. 

4 . Choice  of  Inputs  and  Instrumentation 

Considerable  work  has  been  done  on  the  choice  of  con- 
trol inputs  to  improve  identification  accuracy  of  parameter 
estimates  (see  (iupta-llall  |6]  and  Mehra-Gupta  17|).  Differ- 
ent inputs  provide  different  estimation  errors  for  the  para- 

i o 


4 


meters  of  interest.  Equally  important  is  the  choice  of 
instruments  and  the  location  of  the  chosen  instruments.  In 
many  instances,  however,  the  set  of  instruments  available 
is  fixed  and  it  is  not  practical  to  relocate  the  instruments, 
as  in  the  present  case.  Then  the  choice  of  the  control 
input  is  the  only  design  element  to  improve  the  accuracy  of 
the  parameters. 
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IV.  RESULTS 


4.1  SUMMARY  OF  RESULTS 

The  instrumental  variables  and  the  maximum  likelihood  estima- 
tion methods  are  used  to  identify  the  stability  and  control  coeffi- 
cients of  the  T-2  aircraft.  In  addition,  the  biases  and  noise 
covariances  in  the  instruments  are  also  determined.  One  important 
feature  of  the  computer  program  used  to  process  the  data  is  that  it 
gives  the  covariance  of  the  parameter  estimation  errors  during  the 
identification. 

Flight  test  runs  for  the  low  speed  approach  configuration 
(flight  condition  1)  and  higher  speed  clean  configuration  (flight 
condition  2)  (see  Table  1)  are  processed.  State  noise  is  identi- 
fied in  certain  cases.  The  results  from  different  flight  tests 
at  the  same  flight  condition  are  compared.  Finally,  two  prediction 
runs  are  made,  as  explained  in  Section  4.3,  to  verify  the  identified 
models . 

4.2  LOW  SPEED  APPROACH  CONFIGURATION  (FLIGHT  CONDITION  1) 

The  stability  and  control  coefficients  are  determined  at 
flight  condition  1 (Table  1)  using  response  data  from  runs  2,  3, 

5,  and  6.  These  maneuvers  are  the  responses  to  a two  cycle  sine 
waves  (runs  2 and  5),  a low  amplitude  pulse  (run  3),  and  a high 
amplitude  pulse  (run  6). 

4.2.1  Two  Cycle  Sine  Wave  Elevator  Input  (Run  2) 

The  elevator  input  used  to  excite  the  aircraft  modes  is  shown 
in  Figure  5.  Based  on  the  duration  of  the  run  and  the  input,  there 
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Figure  5 Two  Cycle  Sine  Wave  Elevator  Input  for  Run 


is  only  short  period  excitation.  Since  this  maneuver  is  considered 
a base  line  case,  extensive  processing  was  performed  to  isolate 
various  error  sources  in  the  data.  The  following  paragraphs  detail 
the  effects  on  identified  parameters  for:  (1)  erroneous  acceleration 
measurements,  (2)  non-uniqueness  of  estimates  due  to  identif iability 
problems,  (3)  gust  disturbances,  (4)  noise  in  the  control  measure- 
ment, and  (S)  bias  errors  due  to  out-of-trim. 

Effect  of  Erroneous  Acceleration  Measurement 


The  aerodvnamic  derivatives  are  evaluated  using  an  output 
error  program  in  which  all  parameters  are  assumed  identifiable.  In 
this  run,  measurements  of  normal  acceleration  , pitch  acceleration, 
and  fore-aft  acceleration  are  also  used.  The  estimated  parameters  are 
shown  in  Table  4.  The  actual  response  time  histories  and  computed 
response  time  histories  for  the  estimated  parameters  are  shown  in 
Figure  6.  It  is  clear  that  the  sign  of  fore-aft  acceleration  is 
incorrect  in  the  data.  Another  run  was  made  with  changed  sign  of 
the  fore-aft  acceleration.  The  parameter  estimates  are  shown  in 
Table  4 and  corresponding  time  history  matches  in  Figure  7.  The 
sign  problem  is  solved.  Notice  that  the  short  period  parameters  are 
not  very  different  even  when  an  incorrect  sign  is  taken  for  fore-aft 
acceleration.  Table  5 shows  the  bias  and  white  noise  errors  in  the 
instruments . 

Rank  Deficient  Solution  to  Avoid  Identifiability  Problems 

A rank  deficient  solution  is  used  in  which  all  but  the  short 
period  parameters  are  discarded  since  there  is  insignificant  phu- 
goid  excitation.  Measurements  of  angle-of-attack,  pitch  rate, 
pitch  angle  and  normal  acceleration  are  used  for  identification. 

Table  6 shows  the  wind  tunnel  values,  the  start  up  values  for  maxi- 
mum likelihood  and  the  final  identified  values  assuming  no  process 
noise.  The  corresponding  plots  of  actual  and  predicted  time 
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Table  4 


T-2  Parameters  Using  Output  Error  and  Assuming 
All  Parameters  Identifiable 


corrected  SIGN  POR 
PORI: -APT  ACCELERATION 


INCORRECT  SIGN  POR 
PORE- APT  ACCELERATION 

PARAMETER 

ESTIMATE: 

STANDARD 

DEVIATION 

2 

-0.994 

0.048 

a 

z 

-0.00135 

0.0013 

u 

Z 

1.388 

0.39 

q 

x 

8.76 

3.62 

a 

X 

-0.0168 

0.00689 

u 

M 

-4.42 

0.134 

a 

M 

0.00373 

0.00344 

u 

M 

• 

»— * 
oo 

0.068 

q 

Z«e 

-0.03 

0.1 

X«e 

0.3 

3.46 

M6e 

-8.42 

0.131 

K 

1.0 

a 

(fixed) 

K l 

0 

a 

(fixed) 

ESTIMATE 


-0.965 

-0.0012 

1.23 

-15.82 

0.0329 

-4.52 

0.00434 

-1.863 

-0.07 

-0.3 

-8.562 

1.0 

(fixed) 

0 

(fixed) 


(All  values  in  units  of  ft,  rad,  sec) 

Table  5 

Measurement  Noise  Statistics 


STANDARD 

DEVIATION 


0 .03 
2.52 
0.00427 
0.13 
0 . L00 27 
0.059 
0.085 
1.99 
0.121 


MEASUREMENT 


a 

u 


UNITS 


rad 

ft  sec’ 
rad  sec 
rad 

ft  sec 
rad  sec 


BIAS 

RMS  OF 
WHITE  NOISE 

-0.00301 

0.00642 

-2.41 

3.68 

-0.00113 

0.00535 

0.00933 

0.0117 

-1.19 

0.748 

-0.00102 

0.0636 

Figure  6 Time  History  Plots  for  Run  2 with  Incorrect  Sign  for 


Table  6 

Identified  Derivatives  Using  Rank  Deficient  Solution  - Run  2 

(Approach  Configuration) 


WIND 

TUNNFl. 

START- III*  VAl.lll: S 

PAHAMFTFR 

1-OK  MAXIMUM 
1,1  Ki  1. 1 HOOD 

OUTPUT  FUROR 

PROCFSS  NO  I Sli 

Za 
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-1.046 

-0.903  v 0.032 

-0.873  ♦ 0.0266 

Zu 

0.0 

Zq 

1.0 

1.0 

1.0 

y 

0.0 

Xa 

-23.2 

0.0 

Xu 

-0.046 

0.0 

Xq 

0.0 

0.0 

Ma 

-4.59 

-6.0 

-4.323  ♦ 0.0588 

-4.26  ♦ 0.0922 

Mu 

0.006 

0.0 

-1.42 

-0.48 

-1.658  ♦ 0.103 

-1.543  ♦ 0.114 

! uc 

0.5 

0.5 

Z6e 

-0.102 

0.912 

-0.5724  + 0.787 

0.215  ♦ 0.039 

Xfie 

0.0 

0.0 

M«c 

-9.63 

-8.3978 

-7.445  ♦ 0.196 

-7.32  ♦ 0.33 

Ka 

1.0 

A 

xt /v 

a a 

0.0 

* 

Q 

0.0 

0.000175+0. 0000288 

°a 

0.0063 

0.0091 

0.012S 

°q 

0.0055 

0.0125 

0.0153 

°e 

0.0118 
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0.0170 

°az 

0.71 

4.26 
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(All  values  in  inits  of  ft,  rad,  see) 
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Figure  8 (Concluded)  Ti 


histories  are  in  Figure  8.  Only  five  parameters  are  enough  to 
obtain  a good  fit  to  the  measurements  and,  therefore,  are  the 
only  significant  parameters. 

Process  Noise  Estimation  to  Remove  Gust  Effects 


Next,  the  process  noise  effects  arc  included  in  the  rank 
deficient  solution.  Table  6 compares  the  parameter  estimates 
obtained  from  run  2 response  with  and  without  the  inclusion  of  the 
process  noise  in  the  aircraft  model.  Figure  9 shows  the  time 
history  plots.  There  are  two  major  effects  of  including  the  pro- 
cess noise.  The  parameter  estimates  are  slightly  different  from 
those  obtained  by  the  output  error  techniques  removing,  partly, 
the  bias  in  the  output  error  estimates.  Secondly,  the  one  o 
on  parameter  estimation  error  is  increased,  providing  a more 
realistic  confidence  bound  on  the  errors.  Notice  that 
the  random  noise  in  the  elevator  deflection  measurement  enters 
the  system  as  process  noise,  though  at  a different  point  than  the 
gust  effects.  However,  since  these  two  points  are  not  orthogonal, 
part  of  the  input  noise  is  reflected  in  the  gust  root-mean-square 
value. 


Since  the  identified  power  spectral  density  is  0.000175  and 
the  break  frequency  in  the  Dryden  model  is  0.5  rad  sec  , the 
root -mean- square  gust  is 


^g^rms 


0.000175  = 
2 x 0. 5 


0.0132  rad 


* 0.75  deg 


(7) 


Random  Noise  in  Input 

It  is  clear  from  Figure  5 that  the  measurement  of  the  eleva- 
tor input  is  very  noisy.  A major  part  of  the  error  is  because  of 
noise  in  the  telemetry  link  and  is  essentially  white.  The  input 
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Figure  9 (Concluded)  Time  History  Plots  for  Run  2 with  Process  Noise 


was,  therefore,  smoothed  as  shown  in  Figure  10.  The  SCIDNT  is 
run  using  Figure  10  as  the  elevator  deflection  input  and  the  same 
measurements  as  before.  Table  7 shows  the  identified  derivatives 
with  and  without  the  inclusion  of  process  noise.  Time  history 
plots  of  the  true  and  predicted  measurements  arc  given  in  Figures 
II  and  12  (output  error  and  process  noise  cases,  respectively). 
The  results  are  not  as  good  as  the  output  error  and  the  process 
noise  case  when  the  input  is  not  smoothed.  It  shows  that  an 
optimal  smoothing  as  outlined  in  Appendix  B is  required. 

Trim  Bias 


As  explained  before,  there  is  a bias  in  the  input  because 
the  trim  value  is  not  exactly  known.  In  this  run,  the  input 
bias  is  included  in  the  model  and  is  identified  by  the  maximum 
likelihood  approach  of  Appendix  A.  The  identified  stability  and 
control  derivatives  are  given  in  Table  9 and  the  time  history 
plots  are  shown  in  Figure  13,  for  the  output  error  case.  The 
identified  bias  in  the  input  is  about  0.7  deg.  These 
parameter  estimates  are  more  accurate  because  such  a bias  in  the 
input  is  quite  important.  The  Importance  of  this  bias  is  quite 
evident  in  the  fit  to  the  0 measurement  (compare  with  Figure  8). 

4.2.2  Two  Cycle  S i ii£  Wave  Input (Run  5) 

Run  5 is  at  the  same  flight  condition  as  run  2, and  the  input 
is  the  same,  but  the  data  length  is  shorter.  An  8 see  long  record 
is  used  for  identification  of  aircraft  stability  and  control 
derivatives.  Again,  a rank  deficient  solution  in  the  output 
error  mode  is  used.  Only  the  short  period  parameters  are  identi- 
fiable. The  parameter  values  arc  shown  in  Table  10,  together  with 
measurement  noise  standard  deviation.  The  time  histories  arc 
given  in  Figure  14. 
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Table  7 

Identified  Derivatives  Using  Smoothed  Input  - Run  2 
- 10  Sec  Data  Length 
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(All  values  in  units  of  ft,  see,  rad) 
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Table  8 

bias  and  Initial  Conditions  - Run  2 
Smoothed  Input,  10  Sec 


OUTPUT  liRROR 

PROChSd  NO  I SI: 

BIASHS 

b 

a 

-0.00326 

-0.00230 

bu 

-5.99 

-6.29 

bq 

0. 00853 

-0.0140 

be 

0.00334 

-0.0391 

b3z 

1.77 

1.52 

INITIAL 

CONDITIONS 

a 

0.0 

0.0 

u 

0.0 

0.0 

q 

0.0 

0.0 

e 

0.0 

0.0 

. . 

RUN  2 OUTPUT  ERROR  RUN  2 OUTPUT  ERROR 


RUN  2 OUTPUT  ERROR  RUN  2 OUTPUT  ERROR 


,7 


Figure  11  (Concluded)  Time  History  Plots  for  Run  2 with  Output  Error 

(Smoothed  Input) 
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Figure  12  (Concluded)  Time  History  Plots  for  Run  2 with  Noise 

(Smoothed  Input) 


Table  9 

Identified  Parameters  for  Run  2 
with  Trim  Bias  Removed 


PARAMETER 

WIND 

TUNNEL 

OUTPUT 

ERROR 

Za 

-0.974 

-0.985 

Z 

u 

-0.011 

* 

Z 

q 

1.0 

1.0 

Y 

0.0 

* 

Xa 

-23.2 

* 

Xu 

-0.046 

* 

Xq 

0.0 

* 

Ma 

-4.59 

-5.81 

Mu 

0.006 

* 

% 

-1.42 

-1.52 

U) 

c 

* 

Z«e 

-0.102 

-0.17 

0.0 

* 

M6e 

-9.63 

-7.95 

Ka 

1.0 

*aVV 

0.0 

Q 

* 

0 

a 

0.0097 

0 

q 

J 

0.021 

°o 

0.040 

°az 

4.5 

f A I 1 values  in  units  of  ft,  sec,  rad) 
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Figure  13  Time  History  Plots  for  Run  2 with  Output  Error 

(Trim  Bias  Removed) 
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Figure  13  (Concluded)  Time  History  Plots  for  Run  2 w’th  Output  Error 

(Trim  Bias  Removed) 


Table  10 

Aerodynamic  Derivatives  With  and  Without  Gust  F.ffects 
Doublet  Input  - (Run  5) 


i 

Plight  Condition  1 
V trim  = 255.7  ft.  sec  a trim 

= 5.2° 

WIND 

START  INC 

OUTPUT 

PROCESS  NOISE 

' 

PARAMETER 

TUNNE.i. 

VALUES 

ERROR 

♦ 

i 

Za 

-0.974 

-0.903 

-0.9928  ♦ 0.076 

-0.743  ♦ 0.123 

• 

Zu 

-0.011 

0.0 

1 

1.0 

1.0 

i 

q 

Y 

0.0 

0 0 

Xa 

-23.2 

0.0 

Xu 

-0.046 

Xq 

0.0 

I 

Ma 

-4.59 

I 

-4.706  ♦ 0.0559 

-4.309  ♦ 0.1075 

Mu 

0.006 

0.0 

M 

-1.42 

-1.658 

-1.565  ♦ 0.102 

-1.435  ♦ 0.0604 

q 

u 

0.5 

c 

Z‘e 

-0.102 

-0.572 

-0.354  ♦ 0.121 

-0.407  ♦ 0.1204 

i 

1 

X‘e 

0.0 

0.0 

4 

M«e 

-9.63 

-7.445 

-7.976  ♦ 0.347 

-7.2S  ♦ 0.353 

Xa 

i.O 

WV 

0.0 

Q 

0.0 

1.9  x 10*4  ♦ 2.2  x 10'5 

i 

r 

°a 

0.00911 

0.00468 

0.00412 

o 

0.0125 

0.0112 

0.0093 

q 

°e 

0.0137 

0.00742 

0.005 

°ai 

4.35 

3.34 

3.23 

J 

(All  values  in  units 

of  ft,  sec. 

rad) 
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luded)  Time  History  Plots  for  Run  5 Using  Output  E-.  ror 
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Figure  15  Time  History  Plots  for  Run  5 with  Process  Noise 


ms/u  in 


Figure  15  (Concluded)  Time  History  Plots  for  Run  5 with  Process  Noise 


The  maximum  likelihood  method  is  repeated  with  process  noise. 
The  time  histories  shown  in  Figure  15  have  better  fit.  The  apra- 
meter  estimates  are  given  in  Table  10.  The  computed  one  o values 
are  higher  than  with  output  error  approach.  This  shows  that  the 
output  error  method  tends  to  underestimate  the  parameter  cstima- 
tion  errors.  This  is  an  important  feature  from  the  point  of 
evaluation  of  accuracy  of  estimated  parameter  values.  The  break 
frequency  is  kept  fixed  at  0.5  rad/sec.  Since, 


a = -0.5  a + o) 
g g 


in  the  steady  state 


(8) 


(ag)rms  = / Q/  (l  x T.5)  = /~Q“  (9) 

where  q is  the  power  spectral  density  of  the  process  noise. 
Since  q = 1.9  x 10  4, 


(a  ) = 0.0138  rad  = 0.78° 

g rms 

4.2.3  Pulse  Inputs  (Runs  3 and  6) 


O0) 


The  measurements  of  elevator  deflection  input  and  the  result- 
ing aircraft  response  for  run  3 are  plotted  in  Figure  16.  The 
input  is  a short  duration  pulse  and  there  is  little  excitation 
of  the  aircraft.  It  was  suspected  that  this  data  record  is 
not  going  to  provide  accurate  parameter  estimates.  The  data  is 
run  through  SCIDNT  in  the  output  error  mode.  The  parameter  values 
are  given  in  Table  11  and  the  time  history  plots  in  Figure  16. 

The  standard  deviation  on  these  parameter  estimates  is  higher  than 
for  run  2 or  run  5,  signifying  that  these  parameter  estimates 
are  of  poorer  quality.  Nevertheless,  the  response  time  history 
determined  using  the  parameters  is  close  to  the  measurements. 
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Figure  16  Time  History  Plots  for  Pulse  Input  (Run  3)  With  Output  Error 


RUN  3 OUTPUT  ERROR  RUN  3 OUTPUT  ERROR 


(030)  U13N1 


UJ 

n 


Figure  16  (Concluded)  Time  History  Plots  for  (Run  3)  With  Output  Error 


Table  11 

Flight  Condition  1 - Pulse  Input  (Run  3) 
17.75-23.00  Sec  Data  Length 


1’ARAnrn.k 

WIND 
TUN  Mil. 

START  1 NO 
VAl.lll.S 

OUTPUT  1 UKOK 

1 

a 

-0.974 

-0.974 

-0.788  + 0.1245 

Zu 

-0.011 

-0.011 

* 

z. 

1.0 

1.0 

1 .0 

Y 

0.0 

0.0 

* 

X 

-23.2 

- 25.8 

* 

a 

Xu 

-0.046 

-0.046 

it 

*<. 

0.0 

0.0 

* 

M 

a 

-4.59 

-4.59 

-4.594  + 0.1838 

Mu 

0.0006 

0.0006 

it 

% 

-1.42 

-1.95 

-1.649  + 0.1738 

U) 

• « - 

c 

Z«c 

-0.102 

-0.102 

0.308  ♦ 0.0197 

X«c 

0.0 

0.0 

it 

M6e 

-9.63 

-9.63 

-7.084  + 0.330 

Ka 

1.0 

1.0 

* 

«.VV 

0.0 

0.0 

* 

Q 

— 

— 

— 

°a 

— 

0.007 

0.0086 

°q 

— 

0.015 

0.0142 

°e 

— 

0.015 

0.0078 

°aZ 

— 

0.9 

3.9 

b 

a 

-0.00753 

-0.00753 

bu 

— 

— 

-1.92 

b« 

— 

— 

0.00487 

bo 

— 

— 

0.01 

baz 

.... 

— 

2.0 

(All  valur?  in  units  of  fi  , see,  r.ul ) 
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This  procedure  is  repeated  for  run  6.  The  parameter  values 
are  shown  in  Table  12  and  the  comparison  of  the  time  history  of 
the  measured  response  and  the  estimated  response  is  given  in 
Figure  17. 

4.2.4  Summary  of  Results  and  Model  Verification 


Identification  results  for  flight  condition  1 are  summarized 
in  Table  13.  Also  shown  are  the  aggregate  parameter 
values  by  combining  the  estirates  from  runs  2,  3,  5 and  6 in  the 
inverse  ratio  of  their  covariances.  The  results  from  different 
runs  are  close  to  each  other. 

Next,  the  model  is  verified  by  an  independent  prediction. 

The  parameter  estimates  obtained  from  run  2 with  noise  input  and 
process  noise  consideration  are  using  to  construct  a linear  air- 
craft model  at  the  specified  flight  condition.  This  model  is 
used  with  the  elevator  inputs  for  runs  3 and  5 to  predict  the 
response  of  the  aircraft  to  these  inputs.  These  predicted  re- 
sponses are  compared  to  the  actual  responses  measured  by  instru- 
ments onboard  the  aircraft.  Figure  18  shows  the  comparison  for 
run  3 and  Figure  19  for  run  5.  There  is  a faily  good  match  indi- 
cating that  the  model  is  valid.  Notice  that  the  match  would  be 
better  if  process  noise  is  assumed  in  the  prediction  of  response 
for  runs  3 and  5. 

4.3  HIGH  SPEED,  CLEAN  CONFIGURATION  (FLIGHT  CONDITION  2) 

The  flight  condition  2 is  at  10,000  ft  altitude  and  the  speed 
is  679  ft  sec  The  gear  and  flaps  are  retracted.  Ten  runs 
were  made  at  this  flight  condition.  Out  of  these,  runs  10,  11, 

12,  13,  18  and  19  are  processed  using  SCIDNT. 


Table  12 

Identified  Parameter  for  Pulse  Input  (Run  6) 


PARAMP.TIiR 

WIND 

TIINNl'I. 

STARTING 

VAl.UliS 

OUTPUT  liKROR 

Za 

-0.974 

-0.974 

-0.921  ♦ 0.0454 

Z» 

-0.011 

-0.011 

A 

Z 

q 

1.0 

1.0 

* 

Y 

0.0 

0.0 

* 

X« 

-25.8 

-25.8 

* 

Xu 

-0.046 

-0.046 

A 

Xq 

0.0 

0.0 

A 

Ma 

-4.59 

-4.59 

-4.81  ♦ 0.119 

Mu 

0.006 

0.006 

A 

Nq 

-1.95 

-1.95 

-2.56  ♦ 0.112 

wc 

~ “ 

.... 

~ " 

Z«e 

-0.102 

-0.102 

0.183  ♦ 0.0616 

X«c 

0.0 

0.0 

* 

M<c 

-9.63 

-9.63 

-9.38  ♦ 0.239 

Ka 

— 

1.0 

A 

— 

0.0 

A 

Q 

— 

0.0 

A 

°a 

— 

0.0063 

0.0103 

°q 

— 

0.0055 

0.0095 

°0 

0.0118 

0.0018 

°az 

— 

0.71 

6.2 

Zo 

— 

0.0 

-0.022 

M0 

— 

0.0 

0.0711 

ba 

— 

0.0 

0.0143 

"q 

— 

0.0 

-0.00164 

bo 

— 

0.0 

-0.000567 

b«*z 

.... 

0.0 

- 0 . 000354 

(All  values  in  units  of  ft,  see,  rad) 


Figure  17  Time  History  Plots  for  Pulse  Input  (Run  6) 
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Figure  18  Time  History  Plots  for  Run  3 with  Prediction 
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Figure  18  (Concluded)  Time  History  Plots  for  Run  3 with  Prediction 


After  the  runs  for  flight  condition  1 were  processed,  it 
was  quite  clear  that  it  is  important  to  include  the  possibility 
of  out  of  trim  condition  at  the  beginning  of  the  run.  The  effect 
of  the  out  of  trim  condition  is  accounted  for  by  including  the 
initial  conditions.  An  equivalent  method  is  to  consider  bias 
terms  in  the  equations  of  motion.  In  all  runs  at  flight  condition 
2,  these  equation  biases  are  identified.  The  measurement  biases 
and  mean  square  errors  arc  also  identified. 

Table  14  shows  the  inputs  used  to  excite  the  aircraft  in 
each  of  the  runs.  Runs  10  and  11  use  pulse  inputs  and  run 
12  is  a step  input.  Runs  18  and  13  use  specially  designed  input, 
of  different  amplitudes,  supplied  by  SCI  to  NADC  early  in  the 
program.*  In  Run  19,  random  input  spanning  25  sec  of  data  is 
applied . 

Table  14  shows  the  identified  stability  and  control  deriva- 
tives, equation  and  measurement  biases,  and  the  standard  devia- 
tion of  the  root-mean-square  value  of  the  random  noise  on  the 
measurements.  Also  shown  are  the  standard  deviation  on  parameter 
estimation  errors.  The  output  error  method  is  used  in  each  case. 
Figures  20  to  25  show  the  comparison  of  the  measured  and  esti- 
mated time  history  of  the  measurements  of  angle-of-attack,  for- 
ward speed,  pitch  rate,  pitch  angle  and  vertical  acceleration. 

Note  that  the  fit  between  the  measured  and  the  estimated 
time  history  in  each  case  ;s  good.  The  comparison  of  the  standard 
deviations  on  parameter  estimation  error  shows  that  inputs  of 
runs  18  and  19  are  the  best  and  that  of  run  12  is  the  worst.  The 
input  of  run  12  is  so  small  that  it  cannot  be  distinguished  from 


The  input  was  designed  using  techniques  developed  over  three 
years  ago.  Significantly  better  designs  have  been  developed 
since  then  and  arc  reported  in  References  6 and  7. 
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Table  14 

Identified  Derivatives  from  Different  Runs  at  Flight  Conditi 

Speed  679  ft/sec 
Altitude  10,000  ft 


Table  14  (Concluded) 


« 


i 


> 


i 

1 


r-i 

C 

o 


•o 

c 

o 

u 


JC 

oo 


Uh 


CTl 

</) 

c *-> 

3 <J  4-i 
0£  <D 

w o 

*-»  ^ o 
c +j  o 

<U  44 

t-  o 

O CT>  f“< 
4h 

4-1  vO  0> 
•H  T3 

Q *3  3 

1)  +-» 

E O 

o a 4-« 

u in  T-< 

44  < 

</> 

<U 

> 


rt 
> 
• r4 
»- 
a> 
a 

•a 

o 


c 

0) 

x> 


« 


M 


RUM  10  OUTPUT  ERROR  RUM  10  OUTPUT  ERROR 


1 930 1 IIOIUAJT] 


Figure  20  Pulse  Input  (Run  10) 
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Figure  20  (Concluded)  Pulse  Input  (Run  10) 
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Figure  21  (Continued)  Pulse  Input  (Run  11) 
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Figure  21  (Concluded)  Pulse  Input  (Run  11) 
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Figure  23  (Concluded)  Doublet  Input  (Run  13) 
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Figure  24  Special  Input  ("Run  18) 
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Figure  25  Random  Input  (Run  19) 
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the  measurement  noise  in  the  input.  The  identified  parameters 
from  runs  10,  18  and  19  are  close  to  each  other  confirming  that 
these  inputs  are  "good".  On  the  other  hand,  the  parameter  esti- 
mates from  runs  11,  12  and  13  are  scattered  showing  that  the 
inputs  produce  insufficient  excitation  of  the  aircraft. 

The  equation  bias  terms,  anil  Mq,  are,  in  each  case, 

almost  equivalent  (e.g.,  proportional)  to  the  control  inputs, 

Zg  6e  and  Mg  6e.  This  indicates  that  the  control  input  is  out 

V V 

of  trim.  In  other  words,  the  correct  trim  value  of  the  elevator 
deflection  is  not  subtracted  from  the  time  history  to  determine 
the  input  perturbation  from  trim.  In  run  19,  this  error  is 
about  1°.  Such  large  error  in  trim  value,  if  not  accounted  for, 
would  result  in  a considerable  error  in  estimated  parameters. 

Run  18  is  processed  again  including  the  process  noise  but 
excluding  the  equation  error  terms.  Table  15  shows  the  parameter 
values  and  Figure  26  shows  the  comparison  of  the  measured  and 
estimated  responses.  Since  the  excitation  from  the  deterministic 
input  is  high,  there  does  not  seem  to  be  an  improvement  in  para- 
meter estimates  when  the  process  noise  is  included. 

4.4  IDENTIFICATION  OF  PHUGOID  PARAMETERS 

As  mentioned  earlier,  several  of  the  runs  have  excitation 
of  the  phugoid  motion.  The  maximum  likelihood  computer  program 
SCIDNT  is  used  to  estimate  Z^,  Xu  and  Mu  in  run  10. 
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Table  15 

Flight  Condition  2 - Run  18  (See  Table  1) 
9-16  Sec.  Data  Length 


WIND 

TIINNII. 


-2.4(» 

-o.ooois 


START  INC 
t/Ai.iins 


-2.46 

-0.00018 


I'ROri  SS  NO  I SI: 


-2.682  ♦ 0.759 


36.91 


-0.016 


36.91 

-0.016 


0.0017 


-4.96 


-22.5 

0.0017 

-4.96 


33.43  ♦ 1.06 


-3.858  ♦ 0.343 


-0.1786 


■0.1786 


0.194  ♦ 0.040 


52.12 


52.12 


23.72  ♦ 1.38 


0.0001 

0.007 

0.015 

0.015 


0.0000622  ♦ 
0.0000106 

0.0055 

0.0252 

0.0160 


0.00235 


-1.23 


0.00985 


0.0113 


-3.28 


( A 1 1 v.i  1 nr  s 1 n nn  Us  m ft 
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Figure  26  Time  History  Plots  for  Run  18  with  Noise 


The  phugoid  motion  in  run  10  is  started  by  a pulse  elevator 
input,  which  disturbs  the  airplane  from  trim  conditions.  The 
aircraft  is,  then,  allowed  to  go  through  the  phugoid  oscillation 
without  any  further  input.  A direct  identification  of  the  para- 
meters from  the  recorded  data  gave  erroneous  derivatives  and 
poor  time  history  fits  because  of  excessive  noise  in  the  input. 
When  the  input  is  assumed  zero  for  t > 10  sec,  SC1DNT  gave  the 
time  history  fit  of  Figure  27  and  the  phugoid  parameters  as 
follows 


Z = 0.0009289  ♦ 0.0000394 
u - 

Xu  = -0.0175  ♦ 0.001863 

M = 0.01642  ♦ 0.000424 
u - 

The  identified  phugoid  poles  are  at 
uiph  - -0.0259  ♦ 0 . 1 2 2 j 

It  is  necessary  to  have  better  measurements  of  the  velocity,  input 
and  fore-aft  acceleration  to  estimate  phugoid  parameters. 
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Figure  27  (Continued)  Identification  of  Phugoid  Parameters  (Run  10) 
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Figure  27  (Concluded)  Identification  of  Phugoid  Parameters  (Run  10) 
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V.  SUMMARY  AND  CONCLUSIONS 


S.l  SUMMARY 

A total  of  ten  separate  T-2R  maneuvers  have  been  processed  to 
determine  stability  and  control  derivatives  from  data  sets  supplied 
by  the  Naval  Air  Development  Center.  Of  these  ten,  four  are  low- 
speed  approach  configuration  responses  to  sine  waves,  (two)  and 
pulses  (two)  elevator  inputs.  Six  maneuvers  are  high-speed  clean 
configuration  responses  to  small  and  high  amplitude  pulses  and 
step  a specially  designed  input,  and  a random  input. 

Results  for  the  low-speed  maneuvers  are  summarized  in  Table  13 

of  Chapter  IV  and  in  Figure  28.  Identified  values  of  Za,  Ma,  M^, 

Zfie  and  Mge  are  shown  in  Figure  28  for  the  four  runs  2,  3,  5 and 

6 (ref.  input  schematics  at  bottom  of  figure).  Also  shown  are  the 

lo  standard  deviation  estimates  for  each  identified  run  (e.g.,  the 

larger  the  deviation,  the  less  confidence  in  the  estimate.  On  the 

basis  of  the  identified  values  and  their  respective  confidence 
2 

bounds,  a x statistical  average  of  the  estimates  is  indicated  as 
the  overall  best  estimate  (<). 

Classifying  the  inputs  on  the  basis  of  the  overall  accuracy 
and  consistency  of  the  estimates,  the  following  ranking  from  best 
to  worst  is  concluded: 

(1)  Two  Cycle  Sine  Wave  (Run  2) 

(2)  Two  Cycle  Sine  Wave  (Run  5) 

(3)  Pulse  (run  6) 

(4)  Pulse  (run  3) 


Figure  28  Comparison  of  Stability  and  Control  Derivatives 
at  Low  Speed  Approach  Configuration 
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The  variation  of  Z^e  is  notable  over  the  four  runs,  a result  due  to 
the  difficulty  of  identifying  this  vertical  lift  from  elevator  deflec- 
tion coefficient  without  a sudden  (e.g. , high  frequency)  elevator  input. 


Figure  29  summarizes  the  estimated  values  of  stability  and 
control  derivatives  for  the  six  runs  at  high  speed  and  clean  con- 
figuration. The  estimates  for  the  high  amplitude  pulse  (run  10), 
sepcially  designed  input  (run  i8)  and  random  input  (run  19)  are  the 
best.  The  small  amplitude  pulse  (run  11)  step  (Run  12)  and  specially 
designed  input  (Run  13)  are  the  poorest.  Again,  the  estimates  of 
Zge  are  poor. 

Most  of  the  information  regarding  the  true  value  of  Zfie  is 
contained  in  the  input  measurement  and  the  vertical  acceleration 
measurement.  It  is  concluded  that  both  of  these  measurements 
should  have  smaller  instrument  errors  to  improve  the  estimates  of 
Z^e  with  the  inputs  which  were  tested.  Improved  estimation  of  Z^e 
is  indicated  by  the  use  of  specially  designed  inputs  (e.g. , a high 
frequency  elevator  component). 


An  attempt  is  made  to  identify  the  phugoid  parameters  (run 
10).  It  was  necessary  to  assume  that  the  input  is  zero  beyond 
the  short  period  excitation.  A good  time  history  match  was  obtained 
but  the  parameter  estimates  are  only  fair.  Again,  the  determination 
of  the  phugoid  parameters  requires  a better  measurement  of  the 
input  and  the  aircraft  velocity. 

In  summary,  the  flight  test  data  provides  reliable  and  use- 
ful stability  and  control  derivatives  (and  instrumentation  errors 
and  calibration)  if  the  inputs  are  properly  designed  and  the  errors 
in  important  instruments  are  not  large. 
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Figure  29  Comparison  of  Stability  and  Control  Derivatives  at 
High  Speed  Clean  Configuration 
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S.2  CONCLUSIONS 


On  the  basis  of  processing  of  the  response  data  of  highly 
instrumented  aircraft,  the  following  principal  conclusions  are 
made: 


1.  Identification  Algorithm  Requirements:  In  order  to 

achieve  the  desired  accuracy  of  parameter  estimates,  the 
following  capabilities  are  desirable  for  the  identifica- 
tion algorithm: 

(a)  Trim  state  estimation:  It  is  necessary  to  estimate 

steady  state  errors  due  to  out-of-trim  controls  to 
reduce  parameter  estimate  error. 

(b)  State  noise  estimation:  Significant  reduction  in 

parameter  estimate  errors  is  achievable  by  identi- 
fying the  state  noise  (e.g.,  process  noise)  effects. 
In  particular,  it  is  necessary  to  identify  the 
random  error  in  measurement  of  the  control  inputs. 

(In  these  cases,  identification  of  noise  in  input 
measurement  leads  to  mo^e  consistent  parameter 
estimates . ) 

(c)  Isolation  of  Identifiability  Deficiencies:  It  is 
necessary  to  include  a quantitative  method  of  test- 
ing the  amount  of  information  in  a given  data 
record  to  ascertain  if  particular  parameters  can 
be  identified.  To  not  include  such  a method  can 
lead  to  erroneous  conclusions  about  the  "true" 
values  of  nonidentifiable  parameters. 

2.  Input  Design  and  Measurement  Requirements:  Where  para- 

meters are  not  easily  identified  from  pulses  and  doublets, 
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specially  designed  inputs  (for  desired  short  data  lengths) 
or  random  incuts  (where  arbitrarily  long  data  lengths  are 
practical)  offer  significant  improvements  in  parameter 
estimates.  Particularly  difficult  to  identify  parameters, 
such  as  Zge,  may  require  sudden  step  inputs  to  excite 
high  levels  of  initial  transient  responses.  In  some 
cases,  more  accurate  instrumentation  may  reduce  but  not 
eliminate  requirements  for  such  specialized  inputs. 

3.  Verification  of  Identification  Accuracy:  Evaluation  of 

the  accuracy  of  identified  parameters  should  not  be 
based  on  time  history  matches  between  actual  and  identi- 
fied model  responses  alone.  Estimated  confidence  bounds 
should  also  be  used.  Most  importantly,  the  prediction 
test  is  a very  strong  test  of  estimate  accuracy. 

A central  question  which  is  raised  by  these  results  is  the 
relationship  between  the  wind  tunnel  and  flight  data  identified 
estimates.  The  results  seen  in  Figures  28  and  29  demonstrated 
that  wind  tunnel  static  stability  derivatives  (ZQ,  Ma)  are  verified 
by  the  flight  data  (ref.  statistical  average  of  identified 
derivatives  denoted  by  <) . There  is,  however,  less  validation 
of  wind  tunnel  results  for  the  dynamic  and  control  derivatives 

M6e»  z«e' 

It  cannot  be  concluded  that  one  or  the  other  of  the  results 
is  the  more  valid.  It  may  be  noted,  however,  that  the  wind  tunnel 
values,  placjd  in  a simulation  of  the  aircraft  response,  do  not 
match  the  recorded  aircraft  response,  while  the  identified 
estimates  do  (this  is  implicit  in  the  identification  algorithm 
cost  function).  A combined  fight  test/wind  tunnel  program,  of 
larger  scope  than  conducted  here,  is  necessary  to  isolate  the 
causes  of  the  apparent  discrepancy. 


*.«  1 


Finally,  it  is  concluded  that  flight  testing  of  well- instrumented 
aircraft,  such  as  the  T-2,  offers  an  effective  procedure  for  better 
definition  of  aircraft  aerodynamics  and  instrumentation  evaluation. 

This  procedure  is  significantly  enhanced  oy  use  of  advanced  parameter 
identification  algorithms,  and  careful  attention  to  control  inputs. 

Such  efforts  yield  the  base  for  even  more  effective  research  into 
particular  critical  flight  regimes  such  as  those  at  high  angle-of- 
attack. 
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APPENDIX  A 

AIRCRAFT  PARAMETER  IDENTIFICATION  USING  THE 
MAXIMUM  LIKELIHOOD  METHOD 


Consider  the  aircraft  equations  of  motion  in  discrete  form 


x(t+l)  * *x(t)  ♦ Gu(t)  + rw(t) 


(1) 


y (t)  - Hx (t ) ♦ v(t) 


(2) 


where  x(t)  is  n x 1 state  vector,  u(t)  is  p x 1 input  vec- 
tor, w(t)  is  q x 1 vector  of  random  forcing  functions,  y(t) 
is  r x 1 output  vector,  v(t)  is  r x 1 vector  of  output 
errors  and 


E{w (t ) } = 0,  E{w(t)  w^Ct)}  - Q6t  t 

L > i 


E{w(t ) v (t  ) } * 0 


T, 


E{v(t)  } = 0,  E{v(t)  v (x ) } - R6t 

t t T 


where  6 is  the  Kronecker  Delta  function. 


(3) 

(4) 

(5) 


It  is  assumed  that  the  structure  of  the  model  is  known. 

The  vector  of  unknown  parameters  from  x(0),  <+» , G,  r,  H,  Q, 
and  R is  denoted  by  G.  The  problem  is  to  identify  G using  the 
maximum  likelihood  method. 

For  an  aircraft  in  steady  flight,  and  if  the  wind  gust 
statistics  remain  constant,  the  negative  log* 1 ikel ihood  function 
is  given  by 


!>7 


(6) 


J * 4 l v 1 (i)  B’av (i)  ♦ log  I B I ] 
1 i = l 


where  v(i)  are  the  innovations  from  a Kalman  filter  and  B is 
the  innovations  covariar^e  and  are  given  by  the  following  equa- 
tions 


x(k+l)  = <|>x(k)  + Gu(k)  + <f>Kv(k)  (7) 

v(k)  - y (k)  - Hx (k)  (8) 


K - PHTB-1 

B * HPHT  + R (9) 

P « <p  (P-KBKT)<|>T  + TQrT 

In  the  computer  implementation,  the  negative  log- likelihood 
function  is  minimized  using  a Gauss-Newton  gradient  procedure. 
Minimization  of  (6)  with  respect  to  elements  in  matrix  B gives 
(unknowns  in  R can  be  mapped  into  unknowns  in  B) 

i N T 

& - £ l v(i)  v‘(i)  (10) 

n i = l 


For  other  parameters  0,  the  first  and  second  gradients  of 
J are,  approximately, 


N T . 

Y ( i ) ft'1 
i = 1 


3 v ( i ) 
30 


3J  _ 
30 


(11) 


N 

E 

i = l 


(12) 


a2j 

30^ 


9vT(i)  R-1  »v(i) 

30  D 30 


The  gradients  o£  v arc  obtained  using  Equations  (7)  and  (8) 


3v(k) 

30i 


H ax(k) 

30i 


(13) 


and 


3x(k+l) 

3ei 


, n_u(k) . |i_  Ku(k) 


♦ * If-  '>(*>  * ♦*  <14) 

t h 3K 

where  0.  is  the  i parameter.  The  steady  state  value  of  -2— 

1 oU^ 

is  used  in  Equation  (14).  This  can  be  computed  using  Equations 
(9)  for  each  parameter.  In  the  present  computer  implementation, 
a continuous  system  representation  is  used  to  compute  steady 
state  filter  and  its  gradient. 

In  the  absence  of  process  noise  (wind  gust  effects),  the 
Kalman  filter  gains  and  its  covariance  are  zero.  Therefore, 
Equation  (14)  simplifies  to, 


= H-i(k)  ♦ ^ 


30i 


3x  (k) 
30  • 


3G_ 

30. 


u(k) 


(15) 


The  second  gradient 
which  has  a special 
If 


3 

matrix  -■ is  called  the  information  matrix, 

30  2 

significance  in  parameter  identification. 


(17) 


Standard  deviation  (0.-0.)  = /T7T 

v 1 i'  li 

i.e.,  the  standard  deviation  of  estimation  error  of  the  i*^1  para- 
meter is  the  square  root  of  the  i**1  term  of  the  diagonal  of  Z. 


APPENDIX  B 

MAXIMUM  LIKELIHOOD  ESTIMATION  WITH  BIAS 
AND  RANDOM  NOISE  IN  INPUT  MEASUREMENT 


The  measurement  of  input  signal  is  usually  contaminated  by 
bias,  random  noise  and  other  errors.  These  errors  can  be 
accounted  for  in  the  likelihood  approach.  In  what  is  to  follow, 
we  assume  that  there  is  no  random  excitation  from  gusts  and  the 
input  has  errors  due  to  bias  and  random  noise  only.  Extension 
to  other  cases  is  straightforward. 

Consider  the  discrete  time  representation  of  a system  in 
which  the  state  x follows  the  equation 

x (k+1)  = <{>x(k)  + Gu  (k)  (1) 

There  are  noisy  measurements  of  linear  combinations  of  states 
and  of  inputs 


y (k)  * Hx(k) 

♦ v(k) 

(2) 

u'(k)  = u(k) 

+ b + vu(k) 

(3) 

where  b is  a constant  bias  and  vu(k)  is  white  random  noise. 
Substituting  (3)  in  (1),  we  have 

x (k+1 ) = <*>x(k)  + Gu'(k)  - Gb  - Gvu(k)  (4) 

Equation  (4)  can  be  used  with  measurements  of  liquation  (2)  for 
identification.  Gb  is  the  equation  bias  since  both  G and  b are 
not  time  functions.  The  last  term  in  Equation  (4)  is  a random 
excitation  term  correspond ing  to  process  noise.  Therefore,  even 
though  the  system  docs  not  have  gust  effects,  the  process  noise 
formulation  is  required 


4>x (k)  + Gu'(k)  - Gb  + 4>Kv(k) 


(5) 


x(k+l)  = 

v(k)  - y (k)  - Hx (k)  (6) 

When  the  Newton-Raphson  technique  is  used,  the  gradients  are 
computed  as  before  except  for  the  bias  in  measurement  of  the 
input 


= <J>(I  - KH)  - G (7) 


If  there  were  no  random  noise  in  the  measurement  of  input 


axlkiii  ax^i  . G 


\ 
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